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Under Contract NAS9-X6104, Applied Research Laboratories, The 
University of Texas at Austin (ARL:UT) , was tasked to develop the soft- 
ware required to correlate object scene radar video data frotn the 
APQ-102 radar system. The compressed pulse video data are sampled, 
digitized, and recorded on a wideband tape recorder in the aircraft 
testbed. 

Since the recording format was designed to be compatible with 
reformatting equipment available at ARLUIT (the digital recording inter- 
face equipment, or DRIE) , it was planned to ship all the recorded data 
to ARL:UT so that the Image formation algorithms could be developed and 
the data processed into high resolution radar images. 

A longer term objective of NASA/JSC was to establish the capability 
for reformatting and processing the raw data tapes in their own facility 
but, in the interim, to obtain support from ARLiUT for processing user 
requested radar data. The contractual effort under Contract NAS9-16208 
was, therefore, to perform two principal tasks. 

1. Reformat and correlate selected radar scene data obtained with 
the APQ-102 radar and digital data recording system (DDES) , up to a maxi- 
mum of five images. 

2. Transfer the requisite information for performing the above 
tasks to NASA’s operating contractor, thus establishing the capability 
to complete all data gathering and image generation and display tasks 
within the JSC organization. 
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At the time the contract was initiated, some delay had been 
experienced in developing the software for use on the DDRS data due, 
principally, to lack of DDRS data with which to test the algor itlmis and 
verify the software. Eventually, NASA’s contractor was able to check 
out and render operational the DDRS, but the first of the wideband tapes 
was not provided until the end of October, one month after the contract 
start date. Eventually, the software was developed, but the second 
objective, educating and assisting NASA in the development of an internal 
capability, was carried on simultaneously. 

As a result, when data tapes were finally being generated for 
program investigators, an initial capability had been achieved at NASA 
to generate the Images. A decision was made not to supply those tapes 
to ARL:UT for correlation, but instead to have ARL:UT concentrate on 
refining the software algorithms and assisting JSC in achieving a full 
in-house data processing capability. 


B. Obj ectlves 

Aside from the canceled objective of actually reducing data for 
five images and supplying the results to NASA, efforts were made to 
achieve a final transfer of the software algoritlims to NASA/ JSC for 
implementation on the VAX computer system obtained for the purpose. Com- 
bined with an array processor, the computer gives NASA's contractor a 
powerful tool for the reformatting and processing operations. 

The theory and procedures for image formation were to be conveyed to 
both NASA and Lockheed personnel in tutorial sessions. It was intended 
that the specific algorithms developed by ARL;UT be employed by Lockheed; 
however, once th>a theory was grasped, several nontrivial efforts were 
mounted independently at NASA to process the data from APQ-102 data 
gathering flights using software developed from algorithms based on the 
theory of synthetic aperture radar. 
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Unfoi'tunatCily, the program and subroutines that evolved produced 
Images of relatively low quality, i.e*, they were poorly registered or 
were blurred, ARhrUT's efforts to analyze the problems were unsuccessful 
^ due to inadequate liaison with both NASA and Lockheed personnel, toward 

the end of the contract. In summary, additional work is needed to 
^ investigate the problems with the processing software developed by 

Lockheed , 

In this report, the differences in the two approaches, that 
developed by ARL:UT for NASA, and that developed by Lockheed, are com- 
pared and conclusion$'/ are drawn with respect to NASA's options for 
achieving maxirouro potential from their digital synthetic aperture sensing 
system. 



■ /z::' : 


PRECEDING PAGE BLANK NOT FILMED 


r 


I 

t 

I 



II. 


TASKS 


A. Data Reduction 

The original contract objectives included the processing of five 
radar Images. Data were never provided to ARLtUT for processing, and 
what processing was done was performed by Lockheed at NASA, using rudi- 
meiatary (and incorrect or incomplete) algorithms. The data tapes that 
were required by the contract were not supplied, and a decision was 
eventually made to relieve ARLtUT of this particular requirement. 

B. Technology Transfer 

Following two briefings by ARL;'(IT personnel for Lockheed/NASA 
personnel, one at NASA and one at ARLrUT, a correct analysis of the side- 
looking array radar (SLAR) digital data processing requirements was 
completed by NASA/Lockheed, with the exception of focus requirements and 
aperture weighting effects. However, Lockheed-generated software pro- 
duced radar images of poor quality, and various improvements were 
suggested by ARL:IIT personnel, which were not, unfortunately, based on a 
close liaison with Lockheed. Rather than implementing the ARL:UT developed 
algorithms, the major thrust seemed to be to attempt to correct the Lockheed 
developed software. Although meetings with Lockheed programmers were 
suggested, to try to comprehend their code, particularly the data handling 
routines, such meetings never occurred. 

This report attempts to analyze potential problems in the Lockheed- 
developed software, deficiencies in the algorithms employed, or error 
sources in the parameters and data. Since no significant problems were 
encountered with data supplied to ARL;UT for development of the ground 
signal processor (GSP), it is assumed that the video data are good, and 
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that only errors In parameters (such as starting range sampllrg delay) 
are potential problems, 

The significant changes that have reportedly been Implemented in 
the initial software used by Lockheed are 

1. application of weighting, or a "window**, function to the 
synthetic array, 

2. obtaining and applying the correct STC (sensitivity time control) 
delay to the transmit time, and 

3. applying a focus algorithm to phase weight the data in the 
synthetic array. 

Suggestions that have not been implemented are: 

1, examine the mechanics of synthetic array overlay to ensure that 
near-identical Doppler filters are positioned at the same azimuth for 
each range sample, 

2. add postprocessing algorithms ] 

(1) to use the power in each sample, rather than signal (square | 

root of power) , 

(2) to perform statistical analysis of the overlaid image 
array to determine mean, maximum, minimum, and standard deviation of 
the data, 

(3) to establish a display noise threshold based on image 
statistics, 

(4) after thresholding, to take the logarithm of the data, and 

make gray shade assignments based on a decibel scale and data dynamic i: 

range, and 

(5) to display data in logarithmic form. 
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X III. COMPARATIVE ANALYSIS OF DOPPLER CLUTTER SPECTRUM 

In analyzing the requirements for reducing the raw video pulses 
recorded by the DDRS, the resolution to be achieved in the alongtrack, 
or azimuth dimension, must be considered. In theory, the so-called 
focused array can achieve a resolution of half the real antenna aperture. 
In practice, this is limited by the motion compensation and by the preci- 
sion of the range to the start of sampling, as well as the focus correc- 
tions themselves. 


In the esse of the APQ-102/DDRS system, the delay to the start of 
sampling may be problematic, since the DDRS start time is delayed by an 
arbitrary value from the STC (sensitivity time control) trigger of the 
APQ-102, according to information supplied to ARL:UT by NASA. If the STC 
trigger time depends on the aircraft altitude (J'.n order to maintain a 
constant minimum swath edge nadir angle of 14.2®, for example, in mode 1) , 
the value for the RMIN range delay is readily known from the NASA Earth 
Resources Data Acquisition System (NERDAS) data for the aircraft altitude. 
Figures 1 and 2 indicate the geometry for the two radar modes available. 

ARL:UT, at one point, was told that the STC trigger time could be 
represented as a linear function of altitude H (in feet) for each mode: 

Mode Algorithm Units 

1 0.00195 H + 7.68 psec 

2 0.00165 H + 69.5 psec 

Based on the mapping altitude of 60,550 ft in mode 1, the slant 
range for a nadir angle of 14.2“ is 62,458 ft or a radar delay of 
126.9 psec (see Fig. 1). This is slightly more than the value of 
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X25.8 ysec used in the ARLjUT GSP program for RMIN, which resulted in a 
satisfactory image. Using the algorithm, values are computed as follows. 




Mode 

Altitude 

.(ft)._.. 


Geometry Delay 
(usee) 

Algorithm Delay 
(usee) 

1 

60,000 

14.2® 

125.8 

124.68 

2 

60,000 

45.4® 

173.7 

168.5 

1 

60,550 

14.2® 

126.9 

125.75 

2 

60,550 

45.4“ 

175.2 

169.4 


Although the algorithm has been determined to be in error, and is 
not now in use, the point is that the correct value of delay from trans- 
mit time to the start of sampling must be known, and significant ei" ors 
(greater than 1 psec) must be eliminated. This will be most apparent j 

when the effect of a range error in the algorithms for Doppler filter width 
(array length) and for filter overlay is examined. j 

I 

To assess the overall effectiveness of the correlations, it is ] 

good to first compare the two approaches in representing the usable 

Doppler spectrum. i 

A. NA SA/Lockheed Representation ■ 

' ■ ■ ■ I • ■ ■■ -I ■ - - ■ ^-n-r-rr --fir.i.i . 

Appendix A is the approach to synthetic aperture radar (SAR) processing j 

developed by NASA/Lockheed. Several comments are in order concerning 
these algorithms. First, considering the value of it should be 

pointed out that, for the DDRS initial sample, an additional increment 
of time or delta is set by the operator, ranging from 2 to 

99 jisec. 

A second comment has to do with the maximum unfocused array length 
(see Section III of Appendix A) . The definition is based on the allow- 
able phase shift due to path length distance from one end of the array 
of pulses to be used to the other. This is arbitrarily taken to be 


90°, ±ir/4 radians, with respect to the center of the array. 



IT 


For an unfocused synthetic array of length the backscattered 

nicix 

coherent pulses trom a point target will add constructively with linear 
phase shift in time. Figure 3 illustrates this situation; if we want 
the phase shift corresponding to d to be less than 45® at each end of 
the array, relative to 0® in the middle, then 




9 /L \ 

and since — j » ^ approximated with the first two terms in 

a binomial expansion as 




R /L \ , 

o / max \ I < A 


2 V 2R 


Where the two sides are equal, 


L = 
max > o 


a well-known result. Recalling that the synthetic aperture beamwidth 
is given by 


2L 


the azimuth resolution at range R^ is 


p = 3R = 
a o 2L 


for the nonfocused case. 
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FIGURE 3 

FOR LINEAR FREQUENCY MODULATION 


ARUjUT 
AE-8M73 
CRQ.GA 
12*1 •81 



m 


118 ft 


In the APQ-102 mode 2, Is 136,402 ft, so that 
(from Eq. (3) above) and the maximum number of pulses would be 236 
(2 pulses per foot) for the unfocused array. At the closest range for 
mode 1, <« 62,875 ft, “ 80 ft, axid the maximum number of pulses 

would be 160, If arrays of pulses greater in length than these are to 
be processed, some form of focusing or phase compensation for the end 
points in the array should be implemented. This process is performed in 
the optical processor used by Goodyear with a conical lens. The conical 
lens applies a phase taper as a function of position across the film 
strip in the azimuth direction, and it is also tapered in the range 
direction of the film to account for range variation of the focus function. 

The digital processor implemented by ARL:UT for NASA converts the 
video samples to vectors by applying a unity valued vector with the 
appropriate phase rotation to each of the samples used in the synthetic 
array. A complex Fourier transform results in a complex Doppler filter 
output, and the magnitude is used to form the images. 

In Fig. 4, 1£ a point targat is at y^. and 

R = + y^ + , (6) 

0 T o o ’ 

the two-way phase shift of an RF signal is 

a(t) . - R^(t) . (7) 

where is the time dependent distance from the point target to the 
phase center of the radar. Then 

R(t) = (y^j-vt)^ -f z^ . (8) 
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FIGURE 4 

APQ-1 02 GEOMETRY 
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Expanding tha function R(t) about the point (using Mcbaurin‘s series 


for the derivatives) gives 


R(t) « R, 


1 

2 + *2 


r: 


2 2 
v^t 


+ • . . 


R / R 




The reflected signal from the point target is 


F(t) *• a(t) cos 


w t 
. 0 


- f «(«] 


« a(t) cos 




Ah. D + Ah ( 

o 


2 ' 


4ir 1 


X 2 


0 tv ,2 . 

„2 j R ^ • 

0 


(10) 


This applies during the period of time the point target is illuminated 
by the real I 
gration time, 


by the real beam, i.e., from y^ to y^, and is called the aperture inte- 


^o “ ^1 


The terms in Eq. (10) correspond to a constant phase shift, a constant 
Doppler frequency (linear phase shift with time), and a linear PM 
(quadratic phase shift with time). This last coefficient is 


m., 


2tt v" 


X R. 


1 - 


R 


( 11 ) 


and it is centered on 


w 


j y V 

•'o 


d X R_ 
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There are, of course, higher order terms , which can be ignored in 
the analysis. Putting these results In spherical coordinates, 


^ 0 

— "• sinO cos<}> 
and 


M - ^ ^ (l~3ln^0cos^|^) 


also, 


( 12 ) 


(13) 


w 


d 


V sinG cos<|' 


(lA) 


The Doppler bandwidth is 

I (l-sln^Qco8^({i) I . (15) 

The Doppler spectrum amplitude is, excluding the RF spectrum, 

A(w) « , for |w| - ^ , (lb) 

and is affected by the real antenna beam shape which sets a(t), ns well ns 
the target radar cross section o. For the SLAR case i}!»0, and 

«■ ivvJf 

provided jy -y, | is less than R 6, where B is the 3 dB beamwidth of the 
real antenna beam. 


Thus it is necessary to take into account the FM rate of the Doppler 
returns, m^, since it depends on the range and the angle 0. 


The phase shift ftom the center of the aperture length to one end 


is 


.. 1 "'d^A 

’ 2 ”d \TI “ ~ 


(18) 


To keep A<{i less than ir/4 as a limit, 


vi 2„ 4. Ae 

X R cos 9 g + ^ 

O 0 


2n A ^ Tr 

AR " m 


(19) 


or 


2X 


"o'(l^vdoao)^ 


(20a) 


and, to the other end, 


2X 


R R 4* AR 
o o 


^T^v cos 9^' 


(20b) 


All 


2X 


-R^(R^“AR) 


cosoj 


(21a) 


AR 


2X 


^T^v cosQ^ 


(21b) 


and, assuming 


R > R AR , 
0 0 


|ar| i 


2XR 


(V 
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The depth of the focus is defined as 2R, 


. » « ^ 

(T^v cosOj Ko~^l ® 

If we express the real beam coverage lyQ“y 3 ^l terms of pulses, we have 

N *2|y -y. in feet, so that, for the APQ-102, 

P I o X 


m2 2. 

N cos 6 
P 


For the case of the AFQ-102 operated in mode 1, 7-bit precision, 512 
range bins recorded, and a start sampling delay of 2 psec, we mig’it 
assume a value for of 79,384 ft, X ** 0.1025 ft. For this case, 
0-40.9“ aiTd 

P a 2 - 2 2 ^ 2 . a 3,5 ft for 256 pulses, 55.6 ft for 64 pulses. 

N 

P 


If we relax our criterion for the allowable phase shift to it/2 at each 
end of the synthetic array, so that 


\ ® ‘ (v cose)' 


.V cosO 


). . 


P a a 

D „2 2, 

N cos 6 
P 


the depth of the field is cut in half. 
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It is apparent that for long arrays of pulses, the depth of field 
is very limited, and focusing is advisable. The APO-102 coherently 
detects the bipolar video levels and records the in-phase components 
on film relative to a neutral gray shade corresponding to zero value of 
the in-phase component of the video. These same values are converted 
by the DDRS into 7-, 4-, or 2-bit levels. Miller encoded, and recorded 
on a wideband tape recorder. 

The data can then be focused by applying a phase shift function, 
which effectively converts each value to a complex number. A complex 
Fourier transform produces a filter spectrum containing positive and 
negative coefficients, centered at PRF/4. The filter magnitudes are 
formed by discarding the negative values, which have no physical signifi- 
cance, being merely artifacts of the focusing operation, 

A second and important consideration in the NASA version oi a 
data azimuth (Doppler) correlator is the use of a weighting function. 

The generalized ambiguity function in the crossrange direction has 
-13.2 dB first sidelobes and a -3 dB main lobe (half-power) of 
0.886/B in the Doppler frequency domain, where B is the length of 
the array in time, or half the number of pulses (feet) divided by tlie 
aircraft velocity. In other words, a uniformly weighted array of data 

t 

points results in a sinx/x transformed pattern. To improve the sidelobe 
structure (reduce it), thus reducing sidelobe interference, aperture 
weighting is employed. Table I'*' gives some typical performance parame- 
ters . 


Use of weighting results in main beam broadening and consequent 
decreased resolution. This broadening must be taken into account when 
performing coherent integration of corresponding spectral filters from 
successive data arrays. Although the NASA/Lockheed correlator uses 
the Kaiser weighting function, it is not clear that the beam broadening 
effect is taken into account in the algorithms in the software that 
perform the overlaying, or Integration. 


19 


TABr.E I 

PERFORMANCE FOR VARIOUS FREQUENCY-WEIGHTING FUNCTIONS 




Weighting function 

Ped- 
estal 
height 
H, % 

SNR 

loss, 

db 

Main- 
lobe 
width, 
— 3 db 

Peak 

sideiobe 

level, 

db 

Far 

sideiobe 

fallofT 

Uniform 

100 

0 

0.886/R 

-13.2 

6 db 'octave 

Doiph-Chobyshev 

< 1 » s • s 

* • • » 

1.2/B 

-40 

No decay 

Taylor (A "> 8) 

n 

1.14 

1.25/fl 

-40 

6 dbf/octave 

Co3ine>sauared plus pedestal: 
H + a - H) cos' (ir//B) 
a, Hamming 

8 

1.34 

1.33/fl 

-42.8 

6 db/octavc 

b, 3: 1 “taper ratio” 

33.3 

0.55 

1.09/fl 

-25.7 

6 db, ‘'octave 

cos* 

0 

1.76 

1.46/R 

-31.7 

18 db/octave 

cos* {w//B) 

0 

2.38 

1.66/fl 

-39.1 

24 db>'octave 

cos* (ir//B) 

0 

2.88 

1.94/R 

-47 

30 db^octave 

Triangular; 1 — 2l/|/B 

0 

1.25 

1.27/R 

-26.4 

12 db. octavo 


t In the region |(1 > 8/fl. 


(This table was reproduced from Radar Handbork . M. I. Skolnlk (Ed,).J 
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In datermlning the azimuth resolution, the Important consideration 
Is the shape of the main lobe of the Doppler filter, determined by 
convolving the ''window” function and the uniform distribution of length 
equal to the observation time. The resulting filter shape using, for 
example, the Taylor weighting function with -30 dB first sldelobes, is 
calculated to broaden the -3 dB response by a factor of 1,27. Figure 5 
Illustrates the situation. Figure Hu) shows the -3 dB resolution of 
the filter shape for an unweighted (uniform) collection of sampled data 
(radar pulses) of a given length. Figure 5(b) shows what happens to 
the spectral response after weighting the samples in the time series. 
Because of the finite length of the sample, an array of these responses 
is produced, as illustrated in Fig, 5(c). The response spacing is set 
by the observation time, or alternatively, by the number of pulses times 
the interpulse period or sampling interval. 

Thus, except for a unique combination of observation time and 
window beam broadening factor, the effective linear measure of real 
estate corresponding to a, called RESA (the filter width at the -3 dB 
points), will not be equal to the Doppler sample spacing. The sampling 
ratio (SR), defined as SR"RESA/ sample spacing, is a measure of the 
information content in the image. Specific values of SR should be 
obtained based on user requirements and/or display criteria. 

This means that, in calculating the coverage of each successive 
array of pulses at a given range, the distance moved by the real antenna, 
equal to half the number of pulses in feet, must be converted to an 
equivalent number of Doppler resolution cells for a particular range. 

Concerning the analysis of the problem by NASA/ Lockheed (Appendix A), 
the following comments are In order. First, it is intended to restrict 
the size of arrays processed so that focusing is not necessary. As a 
matter of convenience, however, standard length arrays of 256 or 512 
pulses are Fourier processed, and defocusing takes place. 
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FIGURES 

DOPPLER TRANSFORM FUNCTIONS 







In Section III of Appendix A, the following statement is made: 

•'Since the FFT is in the frequency domain, thus 
each bin (filter component) isj 

A. * PRF/N 

According to sampling theory, the resolution of a sample space 
linear series in the transform domain is no greater than the reciprocal 
of the observation period and, for real valued series, the highest 
recoverable component is half the sample rate* 

This means that by using zero to fill in the pulse sample series, 
any arbitrary resolution may be obtained in the Doppler filter spacing, 
a feature useful in maintaining a given aspect ratio with respect to 
the range dimension of a pixel in the image display. Conversely, using 
more general Fourier transform algorithms, a series of arbitrary length 
may be chosen to obtain a desired filter resolution in the Doppler 
domain. 

Id should be pointed out that the "usable Doppler" of Appendix A, 
Section III, is defined as the Doppler bandwidth encompassed by the 
3 dB beamwidth of the real antenna, On the following page, the usable 
Doppler is also defined as PRF/4 t 1/2 the usable Doppler, i.e., the 
3 dB beamwidth Doppler is shifted from zero fre<|nency to center it on 
PRF/4, since this is the output of the APQ-102 receiver detection system 

The following statement is made in Section IV: 

"Since a primary consideration is the focused part of 
the beam at anyone (sic) given time, then, the total 
number of filters which represent it are: 

J = D^/FT ." 
f 

The approach taken is to select, from an unfocused long synthetic 
array transform those Doppler filters in the central portion of the 
real beam coverage, and assume that they meet the criterion for depth of 
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field. Tnls is incorrect, since these filters way contain contributions 
from phase shifted scatterers at the edges of the real beam coverage 
for the beginning and end points of the synthetic array. 

The proper approach, of course, is to either focus the array prior 
to transformation, or avritably restrict its length. 

Attempts to adjust the focal point for the different ranges to 
compensate for the lack of depth of field will not be successful, 
using the NASA/Lockheed approach. 

B, ARLiUT Representation 


In the ARL:UT derived GSP routines, the azimuth filter spacing is 
usually required to be the same as the range spacing at the patch center 
range, and the azimuth resolution is the same as the range resolution, 
RESA-RESR, at each range bin. The range resolution algorithm is the 
same as in Appendix A. It is Instructive to show that the analyses are 
the same, using the notation of Appendix A, Section III. 
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From Section IV, 
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the number of filters in the real aperture. 
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B *R 
w pc 
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; M*FT “ B *R , 

w pc ’ 


(28) 


J = 


FT 


D * J.RT = JB -R /M 
f w pc/ 


therefore, 


M^XR = J‘B *R , 

pc, w pc 


(29) 


and 


M = J 


.5 JSZ . j „ ^ /m . 

X ’ ^ B WR 

W V pc 


Substituting for M, 
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since, in effect, the PRF is two pulses per foot of platform 
travel, 

N ‘2V , N,, 

This shows that essentially the same formulation is being used by the 
ARLjUT approach in determining the unfocused achievable resolution. 

As pointed out above, the improper lmplcmoi\tation by Lockheed /NASA 
results in a defocused Image, Also, use of the approach in the Lockheed/ 
NASA software can result in the introduction of additional errors, The 
first thing that is done using the Loekheod/NASA approach is to calculate 
the PRF as the reciprocal of twice the ground speed in ft/aee; whereas 
this is correct for the first few arrays formed, the ground speed of the 
plane may change 5% during the time required for image formation. This 
error, depending upon the other parameters, could resiilt in much, larger 
errors of filter spacing, and thus image resolution would bo degraded 
with overlay. 

The ARL;IIT formulation is given in the technical report ARL-TR-Sl-El , 
"Development of a Ground Signal froccssor for Digital Synthetic Army 
Radar Data," 22 May 1981s other considerations, such as beam broadening 
versus filter spacing, accuracy of filter overlay, and selection of N^ 
to achieve desired pixel dimensions, are discussed in that report, 
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IV. RESUtTS 


The APQ-102 radar images produced by Che Lockheed algorithms are 
defocusod (smeared) and create blooming problems on a GET display. The. 
blooming is due to the data being supplied in a linear format. The 
concluding sentence in Appendix A is 

"The data referred to is the square root of the real and 
imaginary pair squared amplitude of the pixel," 

The proper representation should be the power in Che pixel, and 
preferably the logarittm of the power, with suitable thresholds to remove 
noise and adjust the dynamic range of the data to the dynamic range 
of the display. Postprocessing could include a statistical analysis of 
the filter date to establish the thresholds and the gi’ay shade 
assignments . 

As discussed in Section III, there are three possible major error 
sources. 

1. Incorrect or inaccurate range data (timing delay) . 

2. Uncompensated synthetic arrays exceeding the allowable length 
for the nonfocused case. 

3. Improper registration or overlay of filters (azimuth linos) due 
to incorrect data handling in the overlay process or error propagated in 
Che computation of various parameters, 

Appendix B, produced by Lockheed personnel, indicates an attempt to 
address this last item, although the approach is purely empirical, or a 
"cut and Cry" type of solution. 

The discussions of antenna beam shape arc particularly confusing, 
since the antenna beam shape is known to be stox/x in azimuth, and has 
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in fact been measured, 'Xhe filter shape of the Doppler filters enu bo 
derived once the weighting function is chosen, or it can bo determined 
from measurements. The discussion in Appendix B seems to indicato a lack 
of understanding of the relationships between range and the iso-Doppler 
lines of the imago. 
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SAR PROCESSING 


I. The following method is employed 
determine the (Rp^) range to the 
center. Clearly 


- 


+ R‘ 


min 


where 


= 2.5 Nm or 10 Nm 
mm 

dependent upon a mode setting of 
the radar 


to 

patch 



The parameter R^^^^ is found by adding to R^^^ 
the following parameter 

Nrb Sj C/2 

where 

Npg - the maximum number of range bins used 
Sj ” the sample interval sample rate in seconds 
C - the speed of light 

2 - a factor used to denote both sending and receiving a signal 

The parameter Rp^ is found by changing the division factor to 4. 
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II. The slant range resolution is found by the following manner 


hi * *^GRND 
'S'fne' 1 

^GRND " '^Sl/sinG 



III. Azimuth Filter determination 

The FFT once evaluated yields a frequency image of the real aperture; 
however a characteristic of the system which must be kept in mind is 
that the doppler data is offset. By definition, the maximum distance 
in the beam which may be considered as focused is: 


°f ^pc 

Since the data is encoded in PCM form at a rate of 2 pulses per foot, 
then clearly the number of required pulses is: 



Assume an FFT may be found (some power of 2) which contains data 
which is < 2 0^ and 

A. 

Np = 2 - A 


The usable doppler bandwidth is: 
n = 2V-C«SinB 

^op X where 

- real beam v\?idth (1.3°) 
C - a conversion constant 
X - wavelength 
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The zero doppler is: 


D » PRF/4 
zero 

The usable doppler is: 

D-D + D /n 
u zero op/ 2 

Since the FFT is in the frequency domain, thus each bin (filter component) 
is: 

4^ = PRF/Np 

The indices which point to the usable doppler are: 

Is = Di|(-)/Af 

Ie = 

IV. Filter Overlaying 

From section 2 the total number of filters representing the real 
aperture is: 

H = + 1 

This is representative of X number of feet per bin on the ground as: 

FT - B -R , 
w pc/m 

Since a primary consideration is the focused part of the beam at anyone 
given time, then, the total number of filters which represent it are: 

0 = D^/p^ 
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In order to line up successive FFT's the transform must slide in time by 
J*FT*2 pulses. This Implies that the overlay 1s incremented J bins In 
azimuth for a processed range bln. 





Fp r 


The basic philosophy is that if you slide X feet in pulses (shift), it 
also shifts X feet in the aperture. Once shifted, the aperture region 
in common with each other may be simply added. 

The resultant summation of these data may then be output as an image pixel. 
(The data referred to is the square root of the real and imaginary pair 
squared amplitude of the pixel). 


APPENDIX B 

SAR PROCESSING (LOCKHEED) 
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SAR PROCESSING 

Two basic problems have been Investigated for the SAR processing. First the 
type or window functions Msed was made optional, and are 

1 . Cosi ne 

» 

2. Kaiser 

3. Taylor 

A characteristic of these functions are that the Kaiser window tends to have 
a smaller standard deviation from the mean than either the Taylor or Cosine, 

The Taylor window provides a smaller standard deviation than the Cosine. 

The second problem addressed involved overlaying to reduce focusing problems 
at the near and far range bins. By making a modification to permit, if desired, 
the number of pulses used for overlaying variable then parametric runs have 
been generated. 

I 

0 CASE 1 - Using a constant 1.3® beam width, and, the same number of 

pulses in the near and far range to overlap a control case was generated. 

0 CASE 2 - Using a constant 1.3® beam width, (the angle to the near and 
far range bins was assumed to be 1.3®) the respective width of the real 
antenna beam width was computed and used. In this case 68/97 pulses 
were slide for the near and far range bins. 

0 C'lSE 3 - Using the results from CASE 2, the width of the beam was in- 
creased at the near bin to 1.5® and the far end was decreased to 1.25®. 
This effectively Increased the pulses to slide during overlaying over 
case 1 at the extreme ends. 
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The algorithm used to compute the real aperture isi 
1 « 2r tan f /2 

for small angles tan 0 « 0 
1 » 2r f/g * 

where V « the beam width of the antenna 

As long as the angle 1 is to the patch center, this does not seem to generate 
any problem. However, to yield a constant slice of real estate indicate that 
the near angle is > than the middle > the far angle. However, it seems some- 
what intuitive that the antenna pattern is not really rectangular in shape, 
but probably better modeled as a trapozoid. This would permit the angle to 
the main lobe to be 1.3®, the angle at the near end > 1,3® and for end « 1.3®. 


Analysis of Results 

0 CA3E 1 - generates a reasonable picture, but the near and far bins 
tends to be out of focus, (smeared) 

0 CASE 2 - smears the .near range bins and tends to focus the far end, 

0 CASE 3 - has positive results with the focusing of data at the near 
end and far end. The near end bin is not as good as in CASE 
1, the far end is better. 

The real aperture length at the near range bin is assumed '\-r.jV and at the 
far end where r^,rg represent range at the near and far end of the 
image. It is assumed that v. 
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